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The Coulomb exchange interaction is the driving force for quantum coherence in quantum Hall 
systems. We construct a microscopic Landau-site Hamiltonian for the exchange interaction in bilayer 
quantum Hall ferromagnets, which is characterized by the SU(4) isospin structure. By taking 
a continuous limit, the Hamiltonian gives rise to the SU(4) nonlinear sigma model in the von- 
Neumann-lattice formulation. The ground-state energy is evaluated at filling factors v = 1,2,3,4. 
It is shown at v = 1 that there are 3 independent soft waves, where only one soft wave is responsible 
for the coherent tunneling of electrons between the two layers. It is also shown at v = 1 that there 
are 3 independent skyrmion states apart from the translational degree of freedom. They are CP 3 
skyrmions enjoying the spin-charge entanglement confined within the lowest Landau level. 



I. INTRODUCTION 

Quantum coherence is a new aspect of quantum Hall 
(QH) systems. Electron spins are polarized sponta- 
neously due to the exchange Coulomb interaction rather 
than compulsively by the Zccman effect. Skyrmions 
arise as coherent excitations ffl, which have been observed 
experimentally^, |3|, Q as quasiparticles. Intriguingly, an 
interlayer coherence may develop spontaneously between 
the two layers and lead to Josephson-like phenomena 
in bilayer QH (BLQH) systems Q. Recent experimental 
results |(| on tunnelling current may well be interpreted as 
the dc-Joscphson current H though still controversial^. 
We expect even the SU(4) quantum coherence due to 
the spin and layer degrees of freedom^. The driving 
force of quantum coherence is the interlayer exchange 
interaction[jlO| pT[ . The exchange Coulomb interaction 
has also been argued to create a new phase, the canted 
antifcrromagnet pZ 
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|14|, in the BLQH system at the 

filling factor v = 2. 

In this paper we analyze the exchange Coulomb inter- 
action to explore the SU(4) coherence in BLQH ferro- 
magnets. We are concerned about electrons confined to 
the lowest Landau level, where the electron position is 
solely specified by the guiding center (X, Y) obeying 



(1.1) 



This brings in the noncommutative Woo algebra|ll| as 
the basic symmetry of the QH system. It implies that 
the electron position cannot be localized to a point within 
the lowest Landau level, and hence the system cannot be 
described by local field theory. We construct an effective 
field theory to describe phys ics wh ose scale is larger than 
the magnetic length £b = s/h/eB. 

The effective Hamiltonian governing the SU(2) coher- 
ence has been derived [[l0| [ll] by making a derivative ex- 
pansion of the Coulomb energy of spin or pscudospin tex- 
tures at v = 1, where the spin stiffness J is explicitly 
calculated as 



J = 



1 
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(1.2) 



It is a straightforward but complicated task to general- 
ize the SU(2) scheme to the SU(4) scheme, because the 
SU(4) extension of the Woo algebra is considerably com- 
plicated than the SU(2) extension 0. 

We overcome the problem by employing an alternative 
formulation. Namely, we construct a Landau-site Hamil- 
tonian by expanding the electron field operator in terms 
of the one-body wave functions of electrons confined to 
the lowest Landau level. Then, the exchange Coulomb 
interaction emerges just as in ferromagnets. An effec- 
tive Hamiltonian is derived by taking a continuum limit 
in the von-Neumann-latticc formulation, where we sub- 



stitute the spin stiffness (1.2) for the exchange integral. 
This is a consistent procedure in the context of the SU(2) 
coherence. We generalize this procedure to study the 
SU(4) coherence. In the SU(4)-invariant limit the effec- 
tive Hamiltonian is given by the SU(4) nonlinear sigma 
model 



H 1 
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(1.3) 



where T a (x) is the isospin field normalized as J2 a Ta,T a = 
3/8 at v = 1 . The present approach allows us to analyze 
QH states at any filling factor. For the sake of simplic- 
ity, we discuss only integer QH states, though fractional 
QH states are treated similarly in the framework of the 
composite fermion theory Jl6{. 

It is our main purpose to make a thorough investiga- 
tion of soft waves and skyrmion excitations supported 
by the exchange Hamiltonian, though some of them have 
been known previously^ , 17 . We examine carefully what 
are the dynamical fields in the BLQH system. There arc 
three degenerate soft waves in the SU(4)-invariant limit, 
among which only one soft wave is responsible for the co- 
herent tunneling of electrons between the two layers. The 
soft modes are Goldstone modes associated with sponta- 
neous breakdown of the SU(4) isospin symmetry. Actu- 
ally, the degeneracy is resolved by the Zeeman effect and 
the tunneling interaction. Namely, the SU(4) symmetry 
is broken explicitly but softly by these interactions, and 
Goldstone modes turn into pseudo-Goldstone modes with 
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gaps. It is also shown at v = 1 that there are three in- 
dependent skyrmion states apart from the translational 
degree of freedom. They are CP 3 skyrmions enjoying 
the spin-charge entanglement confined within the lowest 
Landau level. 



II. QUANTUM HALL FERROMAGNETS 

To elucidate quantum coherence we start with mono- 
layer QH systems. Electrons make cyclotron motions 
under perpendicular magnetic field B. The number of 
flux quanta passing through the system is = 
where S is the area and <&d = 2irh/e is the flux quan- 
tum. There are iV$ electron states per one Landau level 
by neglecting the spin degree of freedom, each of which is 
associated with one flux quantum. We call it the Landau 
site. One Landau site occupies the area S/N$, = 2irl 2 B: 
and may accommodate two electrons with up and down 
spins. The filling factor is v — N/N& with N the total 
number of electrons. We are concerned about physics 
taking place in the lowest Landau level. 

The microscopic Hamiltonian is a sum of the Coulomb 
term and the Zeeman term, 

Hc = lJ d 2 xd 2 yV(x - y)p(x)p(y), (2.1) 

H Z = -±A Z J d 2 x [pt ( X ) _ p i (*)] , (2.2) 

where V(x — y) = e 2 /4ire\x—y\ is the Coulomb potential; 
p a {x) = i[>' T '(x)ip' T (x) is the electron density with the 
spin index a =|,|; p(x) = pHx) + p l (x); A% = \g* p B B\ 
is the Zeeman gap with g* the magnetic g-factor and ps 
the Bohr magneton. 

We expand the electron field operator in terms of the 
one-body wave functions ipi(x) in the lowest Landau 
level, 



(2.3) 



where c a (i) is the annihilation operator of the up-spin 
(<r =|) or down-spin (a ={) electron at the Landau site 



{c CT (i),4(j)} = 5ij5 ar , 
{c a (i),c T (j)} = {ct(i),ci(j)} = 0. 



(2.4) 



As is well known fl8, 19 , it is impossible to choose an or- 
thonormal comple te se t of one-body wave functions tpi (x) 
in the expansion ([2.3j). Consequently, the electron field 
ip" (x) does not satisfy a standard canonical anticommu- 
tation relation, {ijj (T (x) , ip a ^ (y)} ^ S(x — y), implying 
that an electron cannot be localized to a point within 
the lowest Landau level. Nevertheless, roughly speaking, 
an electron can be localized into a Landau site, which 
has the area 2-k£ 2 b . Thus, it is reasonable that QH effects 
are described by the Landau-site Hamiltonian. 



Substituting (2.3) into (2.1), we derive the direct and 
exchange Coulomb energies, 



(2.5a) 



H X = 2J2Y1 ^4(040X(iK(j)> (2.5b) 

where n(i) = & a (i)c a (i) is the electron number at site 
i: Uij and Jij arc the direct and exchange integrals, 



U, 



d 2 xd 2 y (p*(x)cp*Jy)V(x - y)<p I (x)<p J {y), 



(2.6a) 



■ij --jdxdy ip*(x)ip*{y)V{x - y)ip i (y)(p j (x). 

(2.6b) 

These integrals are convergent because the wave function 
tpi(x) is "localized" within one Landau site i with area 
2itl 2 B . The sum Y^iUj) runs over all spin pairs (i ^ j) 
just once. 

The spin S(i) is defined at each site i by 

with r a the Pauli matrix. Using the algebraic relation 
£4(i)4(?')cr(*')c«r0') = -2S(i)-S(j) - \n(i)n(j), 

CT,T 

(2.8) 

we rewrite the exchange term as 

Hx = -^JijS{i)-S{j) - Jijn(i)n(j). (2.9) 
(*>i) 

The Hamiltonian has the global 0(3) symmetry: It is 
invariant when all spins arc rotated simultaneo usly . 

At v — 1, in the absence of the Zeeman effect (2J2), the 
spin direction is determined spontaneously to minimize 
the exchange energy. The product S(i)-S(j) takes the 
maximum value, S(i)-S(j) = 1/4, when S(i) = S(j). 
Hence, provided > 0, all spins are spontaneously po- 
larized to minimize the exchange energy, where the direc- 
tion of polarization is arbitrary: S(i) = S for all points 
i but the direction of S is arbitrary. Actually, the direc- 
tion of the polarization is the z axis due to the Zeeman 
effect however small it may be. The exchange interaction 
contributes to the ground-state energy by 

(H x ) g = -2j2 Jij = -JV* E J y = -\ £xN ' ( 2 - 10 ) 

{id) 3 



where 



e x = 2j2A 



(2.11) 
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(a) 




FIG. 1: An electron makes a cyclotron motion occupying an 
area 2n£% and avoiding all others. Spin-polarized electrons 
fill the lowest Landau level at the filling factor v = 1. Their 
configuration is represented by a von Neumann lattice with 
the lattice point identified with the center of the cyclotron 
motion. Lattice points in the nearest neighborhood of the 
point x are designated by open circles numbered by 1. Lattice 
points in the second nearest neighborhood are designated by 
double circles numbered by 2. A square lattice (a) and a 
triangular lattice (b) are examples of von Neumann lattices. 



with i fixed arbitrarily; the sum runs over all sites for 
Jij 7^ 0. It is clear that the loss of the exchange energy 
is ex when one electron is removed from filled Landau 
sites. This is equal to the energy necessary to flip one 
spin in the system (|2.9|). 



At v = 2 we obtain (Hx) 
S(i) = and n(i) = 2 in (| 



-exN® = —\exN since 



III. VON NEUMANN LATTICE 

Because the QH system is robust agai nst density fluc- 
tuations, the direct Coulomb term ( 2.5a ) is irrelevant as 
far as perturbative fluctuations are co nce rned. We wish 
to analyze the exchange inte ractio n (2.9). As we have 
stated, the exchange integral (2.6b) is convergent. Recall 
that we have expanded the electro n fi eld in terms of one- 
body wave functions tfi(x) as in ( |2.3| ). The index i may 
represent the angular momentum in the symmetric gauge 
or the linear momentum in the Landau gauge. When we 
evaluate the exchange integral either in the symmetric 
or Landau gauges, we find a large contribution from a 
spin pair (i, j) even if they are not in the nearest neigh- 
borhood of each other. Furthermore, it is not clear at 
all how the rotational and translational symmetries are 
recovered in these gauges when the continuum limit is 
taken: See Ref. jlTj for instance. 

It is most convenient f or u s to use a set of one-body 
wave functions tpi(x) in (2.3) so that the index i runs 
over a lattice such as a square lattice or a triangular 
lattice [Fig.|jJ with the lattice point being the center of 
the cyclotron motion. We can construct such a lattice 
with the use of the coherent-state representation. 

We adopt the symmetric gauge, where A x = \By and 



hb^b in the lowest Landau level with 

bis 7ir B {x+iY) = T2{ z -£? 



(3.1) 



where z = (x + iy)/2£ B - We introduce an eigenstate of 
the angular-momentum lowering operator 6, 



b\P) = f3\(3)- 

Because b is an annihilation operator, the state 
coherent state by definition, and is given by 



_ Bb f -f3*b 



|0) 



e-W /V 6t |0>, 



(3.2) 
) is a 

(3.3) 



where |0) is the angular-momentum zero state obeying 
b |/3) = 0. The wave function (ppix) = (x\(3) is calculated 



<p (x) 



exp 



1 |2 



V2h 



(3.4) 



where (3 = + It describes an electron localized 
around the point z = (3/ V2- 

The coherent state has the minimum uncertainty sub- 
ject to the Heisenberg uncertainty relation associated 
with the noncommutativity (1.1) between the coordi- 
nates X and Y . The state \(3) corresponds to the classical 
state describing a cyclotron motion around the point 



y 



(3.5) 



^Bx. The angular momentum is given by L 



as follows from ( |3.l| ) and (3.2). Because j3 is an arbitrary 
complex number, an electron may be localized around 
any point. 

We consider the QH state at v — 1 . The system is filled 
up with electrons each of which occupies an area 2ir£^. 
It is reasonable to put electrons on a lattice with the 
unit cell area 27r^. Such a lattice is nothing but a von 
Neumann lattice ]2(| [2ll ^3|, |24|. The states on a von 
Neumann lattice form a minimum complete setpOj, |22] ] 
in the lowest Landau level. Thus, we may expan d th e 
electron field in terms of coherent states <fi(x) as in ([2.3]), 
where i runs over all lattice points. 

We consider a square lattice for simplicity [Fig.|](a)]. 
Lattice points are given by (3 mn = * s pn(m + in) , or 

X m = V2^e B m, Y n = -V2^£ B n, (3.6) 

so that the unit cell area is 2ir£ B . States are given by 

\X m , Y„) = cxp[-7r(?™ 2 + n 2 )] exp[V7r(™ + in)& f ] |0). 

(3.7) 

They are not orthogonal, 

(X m , , Y n , \X m ,Y n ) = exp \-~ [(m - m'f + (n - n 1 ) 2 ] ^ . 

(3.8) 
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The wave function (3.4) reads 

(pmn\X) (x\X mi Y n ) 



1 



exp 



x exp I (ym + xn) J , (3.9) 

which describes an electron localized around the lattice 
point (x,y) = (X m ,Y„). 

IV. CONTINUUM LIMIT 

In the von-Neumann-lattice formulation it is straight- 
forward to t ake the field-theoretical limit of the exchange 
energy (2.9), by letting the lattice spacing a vanish just 
as in a lattice model for ferromagnets. The result- 
ing Hamiltonian describes correctly physical phenomena 
whose typical size is much larger than the spacing a. 

We first analyze the nearest-neighbor terms, for which 
we set Jij = J\. Let the lattice points be specified by 



lattice vectors a a with 
product as 



a a = 0. Wc expand the spin 



x a 



(id) 



where a partial integration was made. 
Hamiltonian (2.9) yields 



(4.1) 

The exchange 



Hf~ jj^afaf) J2diS(x)-diS(x) 



(4.2) 



as the lowest order term in the derivative expansion. The 
ground-state energy is given by (2.1C ) with ex = 2 J\ J2 a - 
We next analyze the second nearest-neighbor terms with 
the lattice vectors b^ 3 , for which we set J,, = Jo- We 



obtain the same formula as (4.2) with the replacement of 
a 13 by . Any lattice points can be treated in the same 
way. 

We explicitly consider a square lattice [FigQa)] as 
a simplest von Neumann lattice, where X/«=i a ? a< j = 
2a 2 Sij and X) a 



0(3) nonlinear sigma model 



J d 2 x. Hence, (4.2) amounts to the 



Hf = 2 J 



d 2 x diS(x)-diS(x), 



(4.3) 



where J = J\ and hx = 8Ji in the ground-state energy 
( 2.10 ) for the nearest-neighbor terms. It is easy to see 
that J = Ji+2 J 2 +4J 3 +- ■ ■ and e x = 8 Ji+8 J 2 +8 J 3 +- • • 
by taking into account all lattice points; the series would 
converge rapidly. When we adopt another lattice such as 



the triangular lattice [Fig.[l](b)] and take the contin uum 
limit, we reproduce the same effective Hamiltonian ( |4.3| ) 
together with the ground-state energy (2.10) but with 
different definitions of J and ex in terms of the exchange 
integrals Jy. 

We determine the pa ram eters J and ex as follows. The 
effective Hamiltonian (4.3) was first proposed to study 
skyrmion excitations where the spin stiffness J was 
identified p5l with 



J = 



1 



16x/2^ 47re£ B 



(4.4) 



The formula has been verified by evaluating explicitly the 
energy of a spin texture Jl^, O. We next estimate the 
parameter ex, by substituting the skyrmion configura- 
tion 



5.) 



r 2 + k 



Sy — - 



S z = 



1 T — K 

2~ 



(4.5) 



into the nonlinear sigma model (4.3). One skyrmion in- 
creases the exchange energy [^J by 



(ffx)sky = 47I-J, 



(4.6) 



which is independent of the skyrmion size k in (^j). In 
its small-size limit (k — ► 0) the skyrmion is reduced to 
a hole |^6| . The resultant system is the QH system from 
which one electron has been removed. It corresponds to 
the loss of the exchange energy ex in ferromagnets when 
one electron is removed. Hence, ex = 47r J. 

Consequently, the effective Hamiltonian is given by 



(4.3), which is appropriate to analyze phenomena whose 
scale is larger than the magnetic length £b ■ The ground- 
state exchange energy is given by 



(H 



x) g 



-2-kJN, 



(4.7) 



together wit h th e spin stiffness (4.4). The effective 
Hamiltonian (4.3) describes the spin wave in the QH fer- 
romagnet at v = 1. The spin wave is a Goldstone mode 
associated with the global 0(3) symmetry spontaneously 
broken. Due to the Zeeman effect the Goldstone mode 
acqu ires a gap and the coherent length is made finite: 
See (fc.29a|). 



BILAYER QUANTUM HALL 
FERROMAGNETS 



We generalize arguments to analyze electrons in the 
lowest Landau level in BLQH systems. The SU(2) pseu- 
dospin structure is introduced by assigning up (down) 
pseudospin to the front (back) layer. One Landau site 
contains four electron states in the lowest Landau level, 
which arc distinguished by the SU(4) isospin index a = 
ff,f|,bt,bf. For instance, a = f j implies that the elec- 
tron is in the front layer and its spin is up. The group 
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SU(4) is generated by the Hermitian, traceless, 4x4 ma- 
trices. There are (4 2 — 1) independent matrices. We take 
a standard basis[E7), A a , a = 1,2, • ■ • , 15, normalized as 
Tr(A a Af,) = 2<5 a b- They are the generalization of the Pauli 
matrices. 

We decompose the microscopic Coulomb interaction 
into two terms, 

H+ = i J d 2 xd 2 yV+(x - y) P (x)p(y), (5.1a) 

H c=\j d 2 xd 2 yV^(x - y)Ap(x)Ap(y), (5.1b) 

where depends on the total density p(x), and Hq on 
the density difference Ap(x) between the front and back 
layers, 



Ap{x) = p^{x)+p^(x) - p^ (x) - p bl {x). 



(5.2) 



The Coulomb term Hq is invariant under the SU(4) 
transformation. 



The electron field ^(x) is expanded as in (2.3) 



A'.I, 



(5.3) 



where c a {i) is the annihilation operator of the electron 
with isospin a at site i. Substituting the expansion ( |5.3| ) 
into the Coulomb term (5.1), we extract the direct and 
exchange Coulomb terms. Because the QH system is 
robust against density fluctuations, the dir ect C oulomb 
term arising from the SU(4)-invariant term ( |5.1a ) is irrel- 
evant as far as perturbative fluctuations are conc erned . 
The direct term from the SU(4)-noninvariant term ( 5.11 ) 
is 



H, 



cap 



iV<e 



(5.4) 



where P, = PJ 



i=l 

Pi at each site and 



£caP ~ W, 



5< 



e d2 / u *{l-erf(d/V2£ B )} 



(5.5) 



with the error function erf(x). Here, P^ (i) is the SU(2) 
pscudospin at site i made of the two component spinor 
(cf|,Cb|) as in (2.7). We call the SU(4)-noninvariant 
Coulomb interaction (|5.1b|) the capacitance term since 
£ cap Pz{i)Pz{i) describes the capacitance energy per one 
Landau site. 

We proceed to study the exchange Coulomb interac- 
tion. For this purpose we define the SU(4) isospin at 
each site i by 



(J c t J J 




(5.6) 



Substituting the expansion ( |5.3| ) into (5.1a), and using 
the algebraic relation 

£4(04C7X(i)c«r(j) = -2T(i)-T(j) \n(i)n(j), 



(5.7) 



we obtain the SU(4)-invariant exchange energy as 



H x = - 4 E 4 ( T d)-T(j) + \n{i)n{j) ] 



(hi) 



1 The exchange integral j£j is defined by ( [2.6b| ) with the 
use of t he Coulomb potential V+(x — y). The Hamil- 
tonian (5.8) takes the minimum value when the prod- 
uct T(i)-T(j) takes the maximum value. It occurs for 
T(i) = T(j), where T{i)-T(j) = 3/8 at v = 1 



It is convenient to decompose the exchange term (5.8) 
in terms of various SU(2) components, 



(i,o) 



H x = ~ 4£ 4 ( S*(i)-S*{j) + S h ( l )-S h (.j 



■HI 



4 E4( pT «- pT w+ pi ( i )- pi w) 



■HI 



(i,0) 



2 E4(X> CT (^ ff c?A 



(5.9) 



where (S\i) ■ S\j)) xy = St(i)Si(j) + S y (i)S y (j), etc., 



and 



P a - (Cf T ,C bT )y J , P a - (C U , C H ) - y^j 

(5.10) 

The e xchange energy due to the SU(4)-noninvariant term 
(|5.1b|) is also evaluated. Combining them wc obtain 

Hx = -4E^ (« f (i)-S f (i) + S b (i)-S b (i)) 



(i,3) V 



2E^(X> ct (*k (■?'))> 



(5.11) 
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where jf- 



The exchange Hamiltonian (5.11) is valid at any integer 



filling factor with common exchange integrals Jy and Jy . 
It is an operator which may act on states possessing vari- 
ous isospins. We may restrict the Hilbert space appropri- 
ately at a specific filling factor. We examine two special 
limits to see that it is reduced to the well-established 
results at v = 1. First, we apply a large bias voltage 
and move all electrons to the front layer. The resulting 
system is dynamically equivalent to the monolayer sys- 
tem with spin. Indeed, by using S z = (n 



,1 . 



)/2 and 

it is easy to s ee that fl5.ll ) is reduced to the 



exchange interaction (2.9) describing the monolayer QH 
ferromagnet. Second, we assume a large Zeeman effect 
so that all spins are forced to be polarized, where the 
system describes the spin-frozen bilayer system. W e now 
use P z = (n l — n h )/2 and n = ?i f + ?i b to rewrite ( 5.11 ) 
as 



-4 J2 JijP,(i)PzU) + 4 (Px(i)Px(j) + Py(i)Py(j)) 
(i,3) V 

- Jijn(i)n(j). 



(5.12) 



By taking the continuum limit as in section [V, the ef- 
fective Hamiltonian is found to be 



Hf 



: 2^ j d 2 X {diP x {x)-diP x (x) + diPy(x)-diPy{x)) 

23 Jd 2 xd l P z {x)-d l P z {x), (5.13) 



where the pseudospin field obeys the normalization con- 
dition P(x) 2 = 1/4 at v = 1; the stiffness J is given by 
(O) and 



J d = \ P d% j d 2 q ^V(q)e-^W cxp[-f q 2 ] , (5.14) 



J 



2 d 



02 
1 B 



* B ( 1 - eri(d/V2£ B ] 



(5.15) 



It agrees with the effective Hamiltonian obtained from 
the Coulomb energy of the pseudospin texture fl(i|, [Tl"|l , 
where J = pa and J d = pe in their notation. 



VI. U(l) GAUGE SYMMETRIES 

The density imbalance ctq = 2(P z (i)) between the two 
layers is controlled by applying a bias voltage. It affects 
the system via the interaction term 



Hi 



is = -eVblaa^ Pz(i)- 



(6.1) 



The tunneling interaction is 

H T = -A SAS J2Px(i), 



(6.2) 



where Z\sas is the tunneling gap between the symmetric 
and antisymmetric states. The total Hamiltonian H to t is 
the sum of the exchange term (5.11) and 



JV<6 



H D = - AzS z (i) + e cap P z (i)P z (i) 



A S AsPx(i) - eV biaB Pz(i) 



(6.3) 



It consists of the Zeeman term, the capacitance term, the 
tunneling term and the bias term. 

Since one electron has four components, we may per- 
form local U(4) transformations to the electron field. 
However, the Hamiltonian is not invariant under most 
of them. The symmetry of the direct interaction Ho is 
a direct product of two U(l) symmetries, U T (l)(g)U^ (1), 



ip^(x) 

4> h ^{x) 



Af3( x ) ( i> il (x) 

\^(x) 



(6.4) 



The exchange interaction Hx breaks this into a single 
U(l) symmetry, 



r(x) 



*(as) 



(6.5) 



This is the exact local symmetry of the total Hamilto- 
nian. It should be emphasized, however, that there is no 
gapless mode because there is no propagating mode asso- 
ciated with it: Indeed, the kinetic term of the would-be 
phase field, ip^ (x)dk4>' T { x ), is absent in the Hamiltonian. 

It is imp ortant to recognize |l7| that the gauge sym- 
metry ( |6.5| ) characterizes the BLQH system. [See (8.14) 
why we call it the gauge symmetry.] To show this, let us 
consider a system where the two layers are separated suf- 
ficiently so that there are no interlayer exchange interac- 
tion {Jf- = 0) nor the tunneling interaction (Z\sas = 0). 
Then, the total Hamiltonian is invariant under two local 
transformations, U f (l) and U b (l), which act on electrons 
on the two layers independently, 



f^(x) 
\ip^(x) 

i^{x) 
i> hl (x) 



\ip ts -{x) 

f^(x) 
{ip h ±(x) 



,i/3(x) 



(6.6) 



We may also consider a case without the interlayer ex- 
change interaction (Jfj = 0) but with the tunneling in- 
teraction (Z\sas 7^ 0). Then, the symmetry fl6.6| ) is bro- 
ken into the symmetry ( |6.5| ). The number of U(l) gauge 
symmetries distinguish various bilayer systems. We come 
back to this observation to examine the dynamical de- 



grees of freedom in Section VIH 
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VII. GROUND-STATE ENERGIES 

We evaluate the ground-state energy. Let us consider 
the case v = 1. Unless (P z (i)) = ±1/2 electrons are not 
localized in one of the two layers but rather expand over 
the two layers. The ground state is the up-spin bonding 
state, which is reduced to the up-spin symmetric state 
in the balanced configuration with (P z (i)) = 0. When 
(P z (i)) = (Jo/2, the exchange Coulomb energy reads 

(H x ) v=1 = -£((1 + ol),h, + (1 - al)4) 



= -2n(J + + <jU-)N<s, 



(7.1) 



where 2 J ± = J ± J d . The ground-state energy is 



--(A z 
- 2n(j + 



AsAs) + ^CrQe ca p — -CoeVbias 

ha 2 J-). (7.2) 



In particular, {H) v=1 =-\(Az + A SAS )N$, - 2kJ+N^ 
in the balanced configuration (ctq = 0). 

We next consider the case v = 2. Since two electrons 
exist in one Landau site, we make the composition of 
(pseudo)spins, \®\ = 0©1. We have two types of states 
within the lowest Landau level; (a) three pseudospin- 
singlet and spin-triplet states (the spin sector), (b) three 
pseudospin-triplet and spin-singlet states (the ppin sec- 
tor). 

The spin sector consists of |f^,b^), -L(|fT,b^) + 

If^b^)), |f^,b^). They are the eigenstates of the total 
Hamiltonian within the sector. The ground state is given 
by |f^,b^), and the ground-state energy is 



-(A z +4irJ)N<s,. 



(7.3) 



The state is stable only in the balanced configuration. 
The ppin sector consists of |f^,f^), ^=(|r,b-'-) — 

|f x ,b T )), |b T ,b x ). Within the sector the total Hamilto- 
nian reads 



-cap 



eVbias — ^SAS 



AW* = -A SAS /V2 









-AsAa/y/2 



-A 



SAS 



2 e c 



eVbia 



(7.4) 



apart from a constant exchange energy. The eigenvalue 
equation is easily solved in the balanced configuration 
with the zero bias voltage (Vbias = 0). The ground state 
is given by 



where tan( 



■cos0|S T ,S x ) +sin6»|A T ,A i 



(7.5) 



,/(2A 



SAS 



4^§as 



e? ap )- Here, 



|S^, S^) and | , A^) are the symmetric and antisymmet- 
ric states. The ground state is no longer a symmetric 
state unless e caD = or d = 0. A certain amount of 



the antisymmetric state is necessarily mixed due to the 
capacitance effect. The ground-state energy is 



4^1as + e 



2 ) 

cap / 



2ir(J + J d cos 2 0) 
(7.6) 



in the balanced configuration. The state is stable also in 
unbalanced configurations. 

The (pscudo)spin composition at v = 3 reads \ ®\ ® 

5 = \ ® \ ® I ' wn ere only the doublet is allowed within 
the lowest Landau level. The ground state is a pseu- 
dospin doublet and spin doublet. It is essentially the 
same as the bilayer state at v = 1. The ground state is 

|g)„ =3 = -L(|fT,fT,b^) + |fT,bT,bi)), (7.7) 



and its energy is 



A$ 



-(A z + A sas )-it(5J+ J d ) (7.8) 



in the balanced configuration. The state is stable also in 
unbalanced configuration. 

At v = 4 all the Landau sites are filled up. The ground 
state is pseudospin-singlet and spin-singlet. The ground 
state is 



|fT,fT, b T,b^), 



and its energy is 



(7.9) 



(7.10) 



The state is stable only in the balanced configuration. 

In the SU(4)-invariant limit, where d — > and J d — > J, 
the exchange energy is reduced to a unified formula 
(H x ) = -2ttJN at v = 1,2,3,4. In the SU(4)- 
noninvariant case we note the following intriguing prop- 
erties. (A) In the "layer basis" , where we take four 
independent one-body states |fT), |fj,), |b|) and |bj,), 
the exchange interaction operates only between the same 
isospin states, i.e., 



(H 



x 



-2nJN 



(7.11) 



between a pair of |ff)s, a pair of |f|)s, a pair of |b|)s and a 
pair of |b|)s: All others vanish. (B) In the "SAS" basis, 
where we take four independent one-body states {S^}, 
\Si), \A\) and |A[), we naively expect that the exchange 
interaction operates only between the same isospin states 
as before, i.e., 



(H x ) = — 7r(J + J d )N 



(7.12) 



between a pair of |S"p)s, a pair of \S])b, a pair of I^T) 8 ; 
and a pair of |^4|)s: Actually there appears also an ex- 
change interaction between different isospin states, i.e., 



{H x ) = -Ti(J-J d )N 



(7.13) 
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between |ST) and \A\), and between \Sl) and \ Ai): All 
others vanish. We recover the naive expectation in the 
SU(4) invariant limit. 

We explain why the exchange term J d does not appear 
in the layer basis but does in the SAS b asis. It arises for 
instance from the term P^ (i)-P^ (j) in ( 5.11 ). We find 



(ft|P T (i)|ft) = (£l|P JL (i)|£l> = (0,0 ) i), 



(bt|P T W|bt) = (bi \P l (i)H) 
while 

(St |P T W|St) 
(At |P T (*)|At) 



-(0,0,-), 



(S| |P^)|Sl> = (-,0,0), 



(7.14) 



(A||P^)|A|) = -(-,0,0). (7.15) 

Because only the x and y components contributes to the 
exchange interaction (5.11), there is no contribution in 
the layer basis but there is in the SAS basis. 



VIII. SU(4) SOFT WAVES 

We investigate the SU(4) soft waves at v = 1, which 
are perturbative excitations supported by the exchange 
interaction. To identify the dynamical degree of free- 
dom we use the composite boson (CB) theory of quan- 
tum Hall ferromagncts |26f| by attaching flux quanta to 
electrons p8| p9| , p0[ . The CB field (jf{x) is defined by 
making a singular phase transformation to the electron 
field ip a (x), 



cf>°(x) 



(8.1) 



where the phase field <d(x) attaches one flux quantum 
to each electron via the relation, EijdidjO(x) = 3>-Dp(x). 
We then introduce the normalized CB field n?{x) by 

^(x) = <j){x)n a {x), (8.2) 

so that the TV-component field n a {x) obeys the constraint 

n\x)n{x) = n^{x)n°{x) = 1. (8.3) 



It follows that 4>{x) = p(x). Because the QH system is 
robust against density fluctuations, as far as perturbative 
fluctuations are concerned, we may set p{x) = po, or 



(8.4) 



where po — N/S is the average electron density. 

We count the number of independent fields. The field 
n(x) consists of 4 complex fie lds, but one real field is 
eliminated by the constraint (8.3). Furthermore, the 
U(l) p hase field is not dynamical due to the gauge sym- 
metry (|6.5[), or 



n a (x) 



ia(x) 



n"{x). 



(8.5) 



See also (S.12). This is the only gauge symmetry in the 
BLQH system. Hence, it contains only 3 independent 
complex fields. Such a field is the CP 3 field ]32|. 

A comment is in order. But for the tunneling interac- 
tion and the interlayer exchange interaction, the symme- 
try group is given by (6.6), or 



n^{x) 
n^(x) 

n h ^(x) 
n°^{x) 



n^(x) 
n ^-(x) 

(n h ^{x) 
\n h ^(x) 



(8.6) 



Because there exist two U(l) gauge symmetries, we have 
a set of two CP 1 fields rather than one CP 3 field. 

We conclude that the dynamical field is the CP 3 field 
in the BLQH system due to the exchange interaction 
and the tunneling interaction. See also Section IX. The 
isospin field T and the CP 3 field n are related by 



T a {x) = n\x)^-n(x) 



(8.7) 



Though there are 15 isospin components, only 6 of them 
are independent. 

Let us first analyze the exchange Hamiltonian in the 
SU(4 )-invariant limit, where the exchange interaction 
(pTq) yields a nonlinear sigma model 

15 . 

Hf = 2J + J2 d 2 x[d k T a (x)} 2 . (8.8) 

a=l J 

The SU(4) isospin field obeys the normalization condition 
T{x) 2 = 3/8 at v= 1. 

Since the independent fields are the CP 3 fields n a (x), 
we rewrite the exchange Hamiltonian (3.8) in terms of 
them. Let us define 



T{x) = Y J T a {x) 



(8.9) 



We then have 



T a0 (x) = -^(<5 a/3 - Nn m {x)n a {x)) (8.10) 

with TV = 4 for the SU(4) isospin field. Using this, it is 
straightforward to derive |5l| from ( pls| ) that 



rreff 



I (fx^djrJ -djn) - (n t -a 7 n)(6| 7 n t -n)}. 

(8.11) 



This Hamiltonian has the U(l) gauge symmetry (|3.5]). 
To see this more explicitly, we rewrite it as (31, 32 



Hf = 2J+ J d 2 x{djn^ + iKjrJy {d 3 n - iKjn), 



with 



Kfj,(x) = -in\x)d ll n(x). 



.12) 
.13) 
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The Hamiltonian ( 8.12| ) is invariant under the gauge tran- 
formation 

n°(x) -» e ia ^n a {x), -» /v M + fl M a(s). (8.14) 

Here, the field is not a dynamical field |J2|, since it is 
an auxiliary field defined by flS.13 ) . 

We study small fluctuations of the SU(4) soft waves in 
the balanced configuration with no bias voltage (Vb; as = 
0). The ground state is an up-spin symmetric state at 
v = 1. It is convenient to use the SAS basis rather than 
the layer basis. The ground state is given by 



(n 



ST n si at „M 



n 



(1,0,0,0). 



.15) 



We expand the CP 3 field up to the first order of fluctua- 
tion fields [Fig.@, 



(n s T,n s i,n A T,n A ^(l,Ci,C2,C3), 



where 



;((Ti(x) + i-di(x)). 



(8.16) 
(8.17) 

(8.18) 



They are canonical fields obeying [see fl8.4| )] 
[&(x),(j(y)} = PQ 1 5 tj 5{x - y). 

It is manifest that p^ai{x) denotes the number density 
excited fro m the ground state |S|) to the i-th. level des- 
ignated by (3.16). The field di(x) is the conjugate phase 
variable. 



We expand the exchange interaction (8.12) up to the 
second order, 



{dk'dif}. 



j+ J 



(8.19) 



This Hamiltonian describes three Goldstone modes as- 
sociated with spontaneous symmetry breakdown of the 
SU(4) isospin symmetry. 

Actually, the SU(4) symmetry is broken explicitly but 
softly by various direct interactions. Important SU(2) 
operators are 



S z 

Pz 

Pr 



#1) 



( <7 2 +i) 2 + a 2 + d 2 ) + 



(8.20) 
(8.21) 
(8.22) 



up to the second order of fluctuation fields. Note that 
(3.21) is an exact formula. Direct interaction terms read 



n ,- EcapPO 2 



2 ■ 



n 2 



(of+^+of + i?!), 



Asaspo 2 
rir = -. (ct 2 



(8.23) 
(8.24) 
(8.25) 



(a) A SAS > A z 
v=4 



v=3 

v=2. 
v=1 











^SAS 









A 4 v=3 



(b) A SAS < A z 
v=4 



b4 



v=2 
V=1 





Ik 


i 






c 2 | 




^SAS 



A i 
B j 

A 4 
B 4 



FIG. 2: The lowest Landau level contains four energy levels 
corresponding to the two layers and the two spin states. They 
are represented (a) for ZisAS > Az and (b) for zAsas < Az- 
The lowest-energy level consists of up-spin bonding states, 
and is filled at v — 1. Small fluctuations are Goldstone modes 
Ci, C2 and C3- 



Taking into account the SU(4)-noninvariant exchange in- 
teraction as well, we find that the effective Hamilto- 
nian is decomposed into three independent modes H = 

^spin + Wppin + ^ipin, where 



J+ 



H spin = —{(d k a 1 ) 2 + (d k $ 1 ) 2 } + 



(8.26) 



and 



■11 

2 



A3ASP0 / 2 , q2\ 



.27) 



and 



J+ 



H ipin = —{(d k <j 3 ) 2 + (d k d 3 ) 2 } 



{A Z + A S As)P0 / 2 , a2\ 
(<T 3 + #3 J 



.28) 



They describe three independent soft waves, whic h arc 
pscu do-Go ldstone modes by acquirin g ga ps. Eqs.( |8.26] ) 



and 



.27 ) agree with the results 

©y| derived previ- 
ously for the spin wave and the pscudospin wave (which 
we call the ppin wave), respectivel y. Th e group SU(4) 
is more than SU(2)®SU(2). Eq. ( |8.28| ) is the Hamil- 
tonian obtained newly for the SU(4) component missed 
in the SU(2)®SU(2) component, which we call the ipin 
mode. It is notable that the exchange interactions for 
the spin and ipin modes are solely determ ined by the 
SU(4)-invariant Coulomb interaction fl5.1a| ). The SU(4)- 
noninvariant Coulomb int erac tion contributes only to the 
interaction Hamiltonian ( 3.27 ) o f the ppin mode. This is 
because the noninvariant term ( 5.1b ) involves only the 
density difference Ap or the ppin mode ai = 2P Z . 

The coherence lengths (correlation lengths) are not in- 
finite because the soft modes are gapful. They are 



£spin — t-B\ 



£ipin = i~B\ 



UirJ+ 



A 7 



4ttJ+ 



^SAS 



(8.29a) 
(8.29b) 
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for the spin and ipin modes. The ground state of the ppin 
mode ( fe.27p is a squeezed state |53), where the coherence 
lengths arc different between the conjugate variables 02 
and $21 



/ 4tt J d 

111 ~~ B y 2\T~ 



Sppi 



^SAS 



4ttJ 



ppin 



^cap + ^SAS 



(8.30a) 
(8.30b) 



It is notable that £p P j n is very large when ZisAS is very 



sppm _ 

small. However, £p pin is quite small since e cap is quite 
large in actual samples. 

It is important that the bias voltage Vt>i as couples only 
with the ppin wave. The ipin wave connects the sym- 
metric state with the antisymmetric state and requires 
the tunneling gap for its excitation, but it is insensible 
to the density difference between the t wo la yers. This 
is a direct consequence of the formula (3.21). Further- 
more, it is easy to check that the electromagnetic field 
couples only with the ppin mode because it does not af- 
fect the spin. Consequently, the pscudospin wave is only 
the one that is responsible to the coherent tunneling in 
the BLQH system. The mode has been argued|| to lead 
to the Josephson effect with charge e. 



IX. CP SKYRMIONS 

Provided the Zeeman effect is small enough, 
charged excitations are skyrmions in monolayer QH 
ferromagnets|l|. They are t opo logical solitons in the 
0(3) nonlinear sigma model (4.3). It should be empha- 
sized that the existence of skyrmions is based on the topo- 
logical reasoning. It is argued as follows. The dynami- 
cal field of the nonlinear sigma model is the 0(3) spin 
field S(x). Since it takes value in the 2-sphere S 2 , the 
topological stability is guaranteed based on the theorem 
^(S 2 ) = Z implying that the second homotopy class of 
S 2 is the set of integers Z = {0, ±1, ±2, • • • }. The the- 
orem is rephrased as ^(CP 1 ) = Z. We now argue that 
skyrmions arise based on the theorem ^(CP 3 ) = Z in 
the BLQH system with the SU(4) coherence. 

We consider a generic excitation in SU(N) QH ferro- 
magnets at v = 1. Here, N = 2 in monolayer QH ferro- 
magnets and N = 4 in BLQH ferromagnets. It can be 
proved H ^6) that any excitation confined to the lo wes t 
Landau level is expressed in terms of the CB field ( |8.2| ) 
as 

<P(x) = VpRn CT (z) = e A ^^(z), (9.1) 

where u){z) is an arbitrary analytic function, and A(x) is 
an auxiliary field obeying 



V 2 A(x) = 2ir(p(x) -p Q ). 



(9.2) 



The holomorphicity of u) a {z) in (9.1) is a consequence of 
the requirement that the excitation is confined within the 
lowest Landau level. 

We solve (jO) for the CP^ 1 field, 



n a (x) = 



Substituting (gT|) and (gj) into (gj) we find 
1 



4?r 



V 2 lnp(x) - p(x) +p Q = J° ky (x), 



where 



J° ky (x) = ±-V 2 hiJ2\" a (z)\* 



(9.3) 



(9.4) 



(9.5) 



With the aid of the Cauchy-Ricmann equation for uj(z) 
in ( |9.3D , this is shown p^j to be the time component of 
the topological current|32] defined by 



1 



= ^ x d v K x {x), 



2n 



with ( B.15 ) . The topological charge is given by 



Qsky = / d 2 xJ' ( 



sky ( X )- 



(9.6) 



(9.7) 



It is conserved trivially, d l _ l J^ ky (x) = 0. 

Eq.(| 9 .3 |) is the generic formula for skyrmions fl32|. 
Eq.(|9.4D implies that the density modulation Sp(x) = 
p(x) — Po is induced around a skyrmion. It follows from 
(P~4|) that 



d 2 x[p{x) - p Q ] 



d 2 xJl y (x) = -Q s 



ky 



1. 



(9.8) 



as implies that one skyrmion removes one electron. 

The key of the topological stability is whether the 
skyrmion configuration ( |9.l| ) can be brought into the 
ground-state configuration by a contin uous deformation 
of the CB field. First, the CP^" 1 field Q with Q sky ^ 
cannot be deformed continuously into the ground-state 
value based on the topological theorem ^(CP^ -1 ) = Z. 
Second, the density p{x) cannot be deformed continu- 
ously into the ground-state value po because in the mid- 
stream of this deformation the field configuration escapes 
the lowest Landau level. Indeed, wc have shown that 



p(x) should obey the soliton equation (9.4) as far as 
it is confined within the lowest Landau level. Conse- 
quently, skyrmions are stable in QH systems because 
^(CP^ -1 ) = Z and the QH system is robust against 
density fluctuations. 



The topological charge (9.7) is determined by the high- 
est power of ui a (z). We find Q s k y = n if u) a {z) — * a a z n 



with 



,cr|2 



7^ 0. The lightest skyrmion has the topo- 



logical charge Q 



sky = 1- It is given by the choice of 
b a with a ■ b — in ( |9.3[ ). The skyrmion 
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field (9.3) transforms under the action of SU(N). Since it 
is specified by two parameters a and b with a b = 0, there 
are N(N—1) skyrmion states apart from the translational 
degree of freedom. If the e nerg y is solel y determined by 
the nonlinear sigma model ( |4.3 |) or (|8.8|), all these states 



are degenerate with the energy given by (4.6). When the 
skyrmion is required to approach a specific ground state 
asymptotically, the parameter a is fixed, and hence there 
are N — 1 degenerate skyrmion states. This is physically 
reasonable since there exists one ground state and N — 1 
excitation states in the lowest Landau level [Fig.|2j. 

Let us review skyrmions in monolayer QH ferromag- 
nets (N = 2). The skyrmion is required to approach 
the spin polarized ground state S = (0,0,1/2) asymp- 
totically, and there is no deg ener acy since N — 1 = 1. 
The skyrmion configuration (4.5) is uniquely give n by 
u) = (z, k) in terms of the CP 1 field. It gives ( |4.5| ) via 
S a = \n) T a n. 

We study skyrmions in BLQH systems (N = 4) 
skyrmion is required to approach the ground state 
asymptotically, and there are 3 degenerate states since 
N — 1 = 3. Typical three skyrmions are given by 



The 
H) 



■'spin 



= (2,«,0,0), 
= (Z,0,K, 0), 



(9.9) 



which we call the spin skyrmion, the ppin skyrmion and 
the ipin skyrmion, respectively. They are essentially 
SU(2) skyrmions embedded in the SU(4) theory. 

The degeneracy of these three types of skyrmions 
is resolved by the Zeeman effect and the tunneling 
interaction. Estimation of their excitation energies 
is straightforward Q and compared with experimental 
data|34|, p5[ |. As is obvious in Fig.|^, it depends on a 
competition between the Zeeman effect and the tunnel- 



ing interaction whether spin skyrmions or ppin skyrmions 
are excited thermally. 



X. DISCUSSION 



We have derived the Landau-site Hamiltonian ( 5.11 ) 
for the exchange interaction in BLQH systems. It is valid 
at any integer filling factor. A field-theoretical Hamilto- 
nian is constructed from it based on the von-Neumann- 
lattice formulation. We may use it to analyze phenomena 
whose scale is larger than the magnetic length £b- We 
have analyzed carefully BLQH states at v = 1. The dy- 
namical field is the CP 3 field because of the U(l) gauge 
symmetry inherent in the system. We have found that 
there are three soft waves and three skyrmions. They 
are excitations from the ground state to three excitation 
levels [FigJD in the lowest Landau level. 

Though there are three types of skyrmions, only the 
lightest skyrmions are excited thermally. They are 
spin skyrmions when the Zeeman gap is small enough 
compared with the tunneling gap, while they are ppin 
skyrmions when the tunneling gap is small enough com- 
pared with the Zeeman gap. 

It is interesting to apply the present results to BLQH 
systems at v = 2. In a forthcoming paper we would 
analyze the predicted canted antiferromagnetic phase |l2[ 
|l3| , |l4| . We would also examine a prediction that one 
skyrmion is composed of two skyrmions H , which seems 
to have some experimental supports 
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